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Abstract 

ON 

| We study at- and out-of-equilibrium dynamics of a single homopoly- 

mer chain at low temperature using molecular dynamics simulations. 
The main quantities of interest are the average root mean square dis- 
placement of the monomers below the theta point, and the structure 
factor, as a function of time. The observation of these quantities show 
i-^j ■ a close resemblance to those measured in structural glasses and sug- 

£h I gest that the polymer chain in its low temperature phase is in a glassy 

phase, with its dynamics dominated by traps. In equilibrium, at low 
temperature, we observe the trapping of the monomers and a slowing 
down of the overall motion of the polymer as well as non-exponential 
^ ' relaxation of the structure factor. In out-of-equilibrium, at low teln- 

et ■ peratures, we compute the two-time quantities and observe breaking 

of ergodicity in a range of waiting times, with the onset of aging. 

The problem of the dynamics of a single polymer chain has been addressed 
in many ways, the simplest theory being the Rouse model 0. When the 
polymer is not a Gaussian chain, but evolves either to a swollen or collapsed 
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state after a change of external conditions, the dynamics can also be studied 
using more or less elaborate scenarios, either phenomenologically (like the 
"sausage" or "necklace" pictures of de Gennes and others 0,0]), numerically 
H] or analytically || |(|. Phenomena similar to polymer collapse can be 
commonly observed as a protein folds into its compact native state structure, 
or unfolds to an expanded coil (though the protein problem is much more 
complex than the homopolymer case); more generally a polymer collapses 
when it is put in a poor solvent, and expands in a good solvent. However, 
testing the theoretical results with experiments has, thus far, proven difficult, 
and one has still little insight as to what happens at a microscopic level during 
such dynamics. It is also expected that the dynamics at low temperatures 
is slowed down by traps created by the entanglements of the chain, and the 
compactness of the structure in the collapsed phase. We address this issue 
here for the first time and show that these effects can be measured through 
the one-time and two-time dependence of the structure factor; we prove that 
aging is an intrinsic feature of the dynamics of a long collapsed polymer 
chain. 

We perform the kinetic study of the collapse of a single homopolymer 
chain after a temperature quench, obtained using molecular dynamics simu- 
lations. We characterize the different steps that occur during such collapse. 
We show how these steps can be interpreted on a microscopic level, especially, 
how to account for the fact that for an intermediate time regime the dynamics 
slows down drastically in a way reminiscent of glassy dynamics. We present 
results for one- and two-time quantities and show that, at low temperatures, 
the system behaves in an out-of-equilibrium way. Finally, we demonstrate the 
similarities of a single homopolymer chain with numerical results obtained on 
polymeric systems, such as polymer melts [0], or glass-formers |16|, [17|], with 
the glassy dynamics observed in a one-dimensional model of homopolymer 
collapse |12|], or in 3D heteropolymer and homopolymer dynamics studied 
via analytical techniques [P~3j]. 

We study a "beads on a string" model of a homopolymer in 3D, consisting 
of N = 256 monomers. We model the monomers as hard spheres of unit mass 
with a square-well attractive potential between every pair of contacts. The 
total potential energy of the polymer is 

i N N-1 

£ =nY, u ij+ J2v ii+ x, (1) 
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is the matrix of interactions between monomers i and j, and 



V; 



0, a < |rj - r i+1 | < a x 

+00, |rj - r i+ i| < a , or \r, - r i+1 | > a x 



(3) 



is the potential energy of covalent bonds, is the radius of the hard 

sphere, and ai/2 is the radius of the attractive sphere (e = 1, a — 8.00 
and ai = 11.28). The #-point temperature for the homopolymer studied 
is Tq ps 1.75. The average radius of gyration of the collapsed globule at 
T = 0.2 < T e is i?c ~ 24.7. 

We employ the discrete molecular dynamics algorithm |TJ| . To control the 
temperature of the homopolymer, we introduce a heat bath of 2746 particles, 
which do not interact with the homopolymer or with each other in any way 
but via collisions [15]. By changing the kinetic energy of heat bath particles 
we control the temperature of the environment. The heat bath particles 
are hard spheres of the same radii as the monomers and have unit mass. 
The temperature units are in e/ks- The time unit (see [|l|] for details) is 
estimated from the shortest time between two consequent collisions in the 
system between any two particles. 

Next, we compute: (A) the root mean square displacement (rmsd) R(t w , t w 
r) of the polymer between the two times t w and t w + r: 
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T7^( r i(^) -ri(t w + r)) 2 
, iV i=i 



(4) 



where ri(t w ) and Yi(t w + r) are the coordinates of the monomers at times t w 
and t w + r. The coordinates at times t w and t w + r are adjusted with respect 
to each other so that: (i) homopolymer's centers of mass overlap at these two 
times; (ii) their relative orientation is modified using Kabsch transformation 
r P~5p - (B) The dynamic structure factor is defined as 

F q {t w , t w + r) = 1/N (p* g (t w ) Pq (t w + r)) , (5) 

where p q {t) = Yh exp(iq • r^t)). We compute this quantity for q = q max 
which corresponds to the maximum of the equilibrium structure factor and 
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is equal in our system to q max = 12/90. Thermal averaging, (. . .), is done 
over 10 initial configurations. 

First, we study the polymer at equilibrium, which is achieved by simulat- 
ing the polymer at fixed target temperature for approximately 10 6 time units. 
In this case, we expect time-translational invariance for R(t w , t w + r) = R{t) 
and F q (t w , t w + r) = F q {r). We simulate the polymer at various temperatures 
ranging from 0.2 to 1.75. We find that there are three distinct regions for the 
-R(t) dependence on r at low temperatures (Fig. [I]). The first region can be 
viewed as a ballistic motion of the monomers until they reach their nearest 
neighbors. Indeed, we find that the slope of R{j) vs r on log-log scale is 
close to 1 in this region. Starting at around r ~ 10, the rmsd dependence 
on r becomes slower. The lower the temperature, the flatter is the plot of 
R(t) vs t in this region. Following the dynamics of a single monomer we 
can explain the slowing down of R(t) by the presence of neighbors which 
trap this monomer in a "cage". To escape the trap, a collective effort of the 
monomer and its neighbors is necessary. Thus, this region can be seen as a 
regime of activated processes that drive monomers beyond their cages. At 
later times the dynamics of the system becomes diffusive, where we expect 
and find R(t) ~ y/r, except for T = 0.2 where equilibrium could not be 
achieved during the simulation time. With the increase of temperature the 
trapping time becomes smaller and we observe a change in the curvature of 
the plot R(t) vs r, where the trapping plateau actually disappears. The 
change of curvature occurs at a temperature around T* ps 0.7. 

In order to compare the dynamic behavior of our homopolymer with that 
of known glassy systems (see e.g. fji| [171), we study the structure factor 



F q (r) at equilibrium for various temperatures (Fig. |2|). The mode coupling 
theory (MCT) for glasses |J, has shown that in such systems, the scaling of 
the structure factor behaves non-exponentially with time and predicts power- 
law or stretched-exponential dependence according to the time regime. We 
find that F q (r) has the properties predicted by MCT at low temperatures (see 
Fig. H): (i) initially F q {r) decays as a power- law function 1 — F q {r)/F q {jS) ~ 
r~ 9 to a "plateau", F q ; (ii) F q {j) approaches the plateau, F q , as a power- 
law function (F q (r) — F q ~ (r/r ) _a ); (Hi) F q {r) departs from the plateau 
as a power-law, F q — F q {r) ~ (r/r*) b (j3- relaxation); (iv) in a larger time 
sale {a- relaxation) F q (r) decays to as a stretched exponential, F q (r) ~ 
exp(-(r/T*r). 
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We find that (i) the exponent g = 1.9 ± 0.2, which is in agreement with 
the expected exponent g = 2 due to the ballistic motion of the particles at the 
time scales equal to the size of a typical cage divided by the typical velocity 
of a particle; (%%) exponent a, 0.3 < a < 0.7, is difficult to determine due to 
uncertainty in F^/F q (0) « 0.8 — 0.9; (Hi) the departure from the plateau, 
Fq k, 0.9, is well fitted by a power-law with the exponent b = 0.3±0.02, which 
is almost independent of T; (iv) the a-relaxation is well fitted by a stretched 
exponential with the exponent (3 = 0.52 ± 0.03, which is also independent of 

T m. 

The sign that a system is of out-of-equilibrium is the violation of time- 
translational invariance of quantities such as R(t w , t w + r) and F q (t w , t w + r), 
called aging [JTIJ. To recover aging phenomena, we study the dependence 
of R(t w ,t w + r) and F q (t w ,t w + r) on t w . We quench the polymer from a 
high temperature (T = Tg) to lower target temperatures T t and then let the 
polymer evolve at T t . We start computing the two-time quantities after some 
waiting time t w after the quench (Figs. [| and f|). 

We find that R(t w , t w + r) and F q (t w , t w + r) do not exhibit dependence 
on t w (see Figs. |3]a and £|a) if we quench the polymer to temperatures above 
T t = Tb pa 0.6. At Tb = Tt ~ 0.6, R(t w ,t w + r) and F q (t w ,t w + r) become 
dependent on indicating a slowing down of the dynamics as t w increases. 
Thus, Tb is the temperature at which we observe the breaking of the time 
translational invariance. At T t = 0.3 (see Figs, and £|b) we observe a 
separation of the R(t w , t w + r) and F q (t w , t w + r) vs r curves with respect to 
waiting time t w . In the large t w and r regime, we show that both R(t w , t w +r) 
and F q (t w , t w + r) can be rescaled as a function of r/t 1 ^ (see Figs |3|b and |]b). 
The exponent \i is equal to 0.45 for T t = 0.3 and is the same for R(t w , t w + r) 
and F q (t w ,t w + r). /i depends strongly on the temperature and it increases 
as the temperature decreases. 

In a compact state (at T t ~ 0.3), the polymer exhibits glassy behavior on 
the time scale accessible to our simulations. The spatial rearrangements of 
monomers are slower due to the existence of energetic barriers. Due to the 
finite length of our homopolymer, we expect that at long times, equilibrium 
is achieved and the time-translational invariance is restored. 

We find that the dynamics of a polymer chain at low temperatures is 
similar to the dynamics of strongly frustrated systems, such as glasses [§], 



and of strongly disordered systems, such as spin-glasses [10|. Equilibrium 



plots similar to Fig. [I] and Fig. [| have been previously found in studies of 
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polymer melts |7[]. Here we also demonstrate that the existence of aging is 
an intrinsic feature of the dynamics of a polymer chain at low temperatures. 
We expect the microscopic mechanisms responsible for the slowing down 
of the dynamics to be similar to those found in glasses, and they can be 
easily identified in our case as geometrical frustration. The two temperatures, 
T* pa 0.7 and T& ~ 0.6, at which trapping starts to dominate the dynamics 
and time translational invariance is broken, are close to each other and may 
well correspond to a unique dynamical transition temperature. The analysis 
of the equilibrium data is consistent with some of the results predicted by 
the MCT theory. From the point of view of aging, we do not know any 
theory for 3D polymer dynamics that would predict ar/tjj, scaling. However, 
this scaling (with /i < 1) can be found in a number of systems, including 



real spin-glasses [[nj, and in models of traps |TT| or in a model of a one- 
dimensional polymer [ fL2|| ; whereas a r/t w scaling is more specific to near- 
mean-field models [|1(], [T3|j . 

A single homopolymer provides a new example of a non-disordered sys- 
tem where geometry alone induces slow relaxation and aging. A number of 
polymeric systems in a compact state (single polymeric molecules, polymer 
melts, collapsed proteins) share the dynamical features described here. Some 
open questions remain, such as the very long time regime of the relaxation 
which is hard to characterize due to long computational times. 

We thank A. Coniglio, W. Kob, S. Mossa, F. Sciortino, and E. I. Shakhnovich 
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Figure 1: The dependence of the rmsd R(t) on r in equilibrium simulations 
at various fixed temperatures for a polymer chain of 256 monomers. Three 
regimes can be observed: (i) ballistic (at small time scales), (ii) dominated 
by traps (at intermediate time scales), and (Hi) diffusive (at long time scales). 
As the temperature increases, the second regime becomes less apparent and, 
at temperatures around T* w 0.7, it disappears. 
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time, x 

Figure 2: The dependence of the structure function, F q {r) on r in equilib- 
rium simulations at various fixed temperatures. We choose q = q max , which 
corresponds to the nearest neighbor shell of the polymer. 
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Figure 3: The dependence of the rmsd R(t w ,t w + r) on r in quenching 
simulations from Tg « 1.75 to (a) T t = 0.6 (no aging observed) and (b) 
T t = 0.3 (aging) for various waiting times, t w . Inset: rescaling of the curves 
for t w = 10 5 , 10 4 , 10 3 , with respect to r/tf^, where \i ~ 0.45. The drastic 
difference between (a) and (b) points towards the existence of a temperature 
at which the polymer dynamics becomes non-ergodic. 
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Figure 4: The same as Fig. § but for the structure factor F q (t w ,t w + r). 
The exponent \i is equal to 0.45. 
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